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Resumen 


CAMPOS-ARANDA, D.F. Estimación de las magnitudes 
asociadas con el rompimiento de presas de tierra O 
enrocamiento a través del método estadístico. Tecnología y 
Ciencias del Agua. Vol. V, núm. 3, mayo-junio de 2014, pp. 
159-172. 


Los estudios de seguridad de presas definen los planes de 
emergencia y las estrategias de mitigación de riesgos en 
sus áreas localizadas aguas abajo. Tales estudios abarcan la 
estimación de las magnitudes originadas por el rompimiento 
de su cortina, las cuales incluyen la geometría de la brecha 
de falla, el tiempo de desarrollo de ésta y el gasto máximo 
de descarga. Existen tres enfoques de estimación de las 
características citadas: (1) el método estadístico, que aplica 
ecuaciones empíricas obtenidas por regresión; (2) el método 
hidráulico, el cual simula numéricamente la formación de 
la brecha de falla, tomando en cuenta los mecanismos 
hidráulicos y de erosión del flujo, y (3) el método experimental 
o de laboratorio que utiliza modelos reducidos o a escala 
real. En este trabajo se aplica el enfoque estadístico basado 
en 28 ecuaciones empíricas que se han establecido como 
las más confiables o de menor incertidumbre. Se describen 
dos aplicaciones numéricas hipotéticas: la primera en una 
presa pequeña que falla por desbordamiento; la segunda en 
una presa de mediana altura, cuya falla se debe a la erosión 
interna que generan las filtraciones; en ambas aplicaciones 
se siguen dos planteamientos de análisis, función de la 
información disponible. Se concluye que la recopilación 
expuesta de 28 ecuaciones empíricas permite establecer los 
valores probables de las magnitudes asociadas con la falla 
de presas con cortinas de tierra o de enrocamiento, así como 


sus intervalos de variación. 


Palabras clave: cortinas de tierra, cortinas de enrocamiento, 
desbordamiento, tubificación, brecha de falla, tiempo de 


formación de la brecha, gasto máximo de descarga. 


Abstract 


CAMPOS-ARANDA, D.E Estimate of the Magnitudes of 
Ruptures of Earth or Rockfill Dams using a Statistical Method. 
Water Technology and Sciences (in Spanish). Vol. V, No. 3, May- 
June, 2014, pp. 159-172. 


Studies of dam safety define emergency plans and strategies to 
mitigate risks in underwater areas. These studies involve the 
estimate of the magnitudes of dam ruptures, which include the 
geometry of the failure breach, the time ¡it took to develop and peak 
discharge flow. Three approaches to estimating these characteristics 
exist: (1) the statistical method which applies empirical equations 
obtained through regression, (2) the hydraulic method which 
numerically simulates the development of the failure breach, taking 
into account hydraulic mechanisms and flow erosion and (3) the 
experimental method or laboratory method which uses models at 
reduced or actual scale. This work applies the statistical approach 
based on 28 empirical equations which have been established as most 
reliable or having the least uncertainty. Two hypothetical numerical 
applications are described— the first for a small dam with failure 
due to overtopping and the second for a dam of medium height 
which failed due to internal erosion that caused filtrations. Both 
applications follow two analytical approaches depending on the 
available information. It is concluded that the collection of the 28 
empirical equations makes it possible to establish the probable values 
for the magnitudes associated with the failure of earth or rockfill 
dams, as well as variation intervals. 


Keywords: Earth dam, rockfill dam, overtopping, piping, failure 
breach, time to develop breach, peak discharge flow. 
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Introducción 


A partir de los años ochenta creció el interés 
en el mundo por la seguridad de las presas, 
teniendo como objetivo fundamental la reduc- 
ción de la pérdida de vidas humanas y de los 
daños materiales aguas abajo causados por su 
falla o rompimiento. En realidad, la falla de 
presas es infrecuente, pero sus consecuencias 
son, en general, catastróficas (NWS, 2002). 

Se han construido más de 800 000 presas 
en el orbe, así como miles de kilómetros de 
diques o bordos de protección a ciudades y 
zonas con potencial económico. La historia de 
su construcción también registra que cientos de 
presas han fallado en los últimos siglos y que 
actualmente cada año se rompen tramos de 
diques o bordos debido a su desbordamiento 
u otras causas (Zagonjolli, 2007). Ejemplos 
dramáticos de esto último son la ciudad de 
Nueva Orleans, Estados Unidos, en agosto de 
2005, bajo la influencia del huracán Katrina, y 
la zona de Chalco en el Valle de México, con la 
falla periódica de los bordos del río o canal de 
La Compañía. 

En general, todas las presas están sujetas 
a envejecimiento y como resultado de ello 
se vuelven más peligrosas, por ser más 
susceptibles de fallar por deficiencias en su 
operación o por desbordamiento asociado con 
su azolvamiento; esto último, principalmente 
en presas pequeñas, que son la mayoría. 
Además, las presas antiguas fueron diseñadas 
con información escasa y construidas quizás 
con criterios menos estrictos que los siguidos 
hoy día. Por ejemplo, las normas de diseño y 
construcción de presas de tierra establecen que 
tanto la cimentación como la cortina cumplan 
con criterios muy precisos de control de 
filtraciones, presiones y estabilidad, así como 
medidas de seguridad hidrológica contra el 
desbordamiento. Sin embargo, la falla de presas 
de tierra no ha disminuido de forma sustancial 
(Singh, 1996). 

Por otra parte, los efectos del cambio climáti- 
co es probable que establezcan condiciones más 


extremas de crecientes en las presas existentes y 
por ello se deben estudiar intentando modelar 
su falla o rompimiento, lo cual involucra dos 
estimaciones fundamentales: (1) el hidrograma 
de egresos, caracterizado por su gasto máximo, 
mismo que es función de la geometría de la 
brecha de falla y del tiempo de desarrollo de 
ésta; y (2) el tránsito de la onda de la creciente 
originada por la falla hacia aguas abajo hasta 
el valle y/o la planicie de inundación (Singh, 
1996; Chinnarasti et al., 2004). 

Para estimar las magnitudes citadas, aso- 
ciadas con la falla de una presa, se tienen tres 
enfoques de aproximación: (1) el estadístico, 
orientado al desarrollo de ecuaciones empíricas 
por regresión; (2) los métodos hidráulicos, los 
cuales simulan de modo numérico la formación 
de la brecha de falla, tomando en cuenta los 
mecanismos hidráulicos y de erosión del flujo; 
y (3) los métodos experimentales o de laboratorio 
que utilizan modelos reducidos o a escala real, 
y que permiten verificar ya sea los resultados 
de los otros dos enfoques o bien generan 
nuevas relaciones empíricas (Chinnarasri et al., 
2004; Zagonjolli, 2007). 

El objetivo de este trabajo consiste en 
exponer las ecuaciones de regresión que se 
han establecido como las más confiables o 
de menor incertidumbre para la estimación 
de la 
brecha de falla, su tiempo de formación y 


de las características geométricas 
el gasto máximo descargado por ésta. Estas 
28 ecuaciones empíricas se integran en tres 
grupos, según su grado de incertidumbre. 
Se describen dos aplicaciones numéricas 
hipotéticas, la primera en una presa pequeña 
que falla por desbordamiento y la segunda en 
una presa de mediana altura, cuya falla se debe 
a la erosión interna que generan las filtraciones 
(tubificación); en ambas aplicaciones se siguen 
dos planteamientos de análisis que son función 
de la información disponible. Con base en 
sus resultados se establecen las magnitudes 
buscadas asociadas con la falla o rompimiento 
de una presa con cortina de materiales térreos, 
incluyendo el intervalo de variación de tales 
aproximaciones. 
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Recopilación de ecuaciones empíricas 
Información sobre fallas de presas 


La cortina de una presa puede fallar de manera 
gradual o súbita. El tipo de falla depende en 
particular del material de la cortina. Por lo 
general, las de concreto fallan al instante, y 
las de tierra o de enrocamiento lo hacen de 
forma gradual. Las cortinas más comunes son 
las de materiales térreos (homogéneas, mixtas 
y de enrocamiento), y sus fallas están más 
documentadas, con duraciones rápidas de 15 
minutos a 3 horas, y con lapsos prolongados de 
3 a 12 horas, y aún más (Singh, 1996). 

Varios investigadores exponen en sus 
trabajos los archivos que emplearon para 
deducir sus ecuaciones empíricas propuestas, 
por ejemplo, Hagen (1982), MacDonald y 
Langridge-Monopolis (1984) Costa (1985, 
1988), Froehlich (1995a, 2008), y Xu y Zhang 
(2009). 


Dimensiones necesarias y estimaciones buscadas 


De manera general, las magnitudes relativas 
al rompimiento o falla de una presa se pueden 
dividir en dos grupos: parámetros geométricos 
e hidrográficos. Los primeros corresponden 
a la brecha de falla, la cual por lo común es 
trapezoidal, con cinco parámetros geométricos: 
ancho máximo (b,,), ancho medio (b), ancho 
del fondo o mínimo (6, ), tirante o altura de la 
brecha (H,), y taludes de sus lados (z). En la 
figura 1 se indican los elementos geométricos 
generales de una brecha de falla. Los paráme- 
tros hidrográficos son el tiempo que dura la 
falla (1) y el gasto máximo descargado (Q,). 
El T, en rompimientos por desbordamiento es 
el lapso desde el inicio del crecimiento rápido 
hacia abajo hasta que cesa la erosión lateral 
de la cortina (Chinnarasri et al., 2004). En las 
fallas por erosión interna debido a filtraciones, 
abarca desde su inicio el colapso de la cresta de 
la cortina y la subsecuente formación lateral de 
la brecha (Xu y Zhang, 2009). 

La aplicación de la mayoría de las ecuacio- 
nes empíricas requiere como datos básicos las 


siguientes dos características geométricas O 
hidráulicas de la presa (NWS, 2002): 


1. Altura hidráulica (H,) o altura de la cortina 
(H.): es la distancia en metros desde el nivel 
del cauce en el sitio de la cortina hasta la 
elevación máxima del agua almacenada. 

2. Volumen de almacenamiento (V.): espacio 
máximo para almacenamiento en mf? en 
la presa bajo el nivel máximo alcanzable 
por el agua, incluyendo cualquier sobre- 
almacenamiento. 


La primera definición establece que H, = H, 
entonces implica tres simplificaciones: (1) que 
el piso o fondo de la brecha de falla alcanza el 
cauce, es decir, que toda la altura de cortina 
se erosiona; (2) que en fallas por tubificación 
el nivel del agua llega hasta la corona de la 
cortina, y (3) que en fallas por desbordamiento 
no se considera la lámina de agua sobre la cresta 
de la cortina. Cuando el fondo de la brecha de 
falla no llega al cauce y se conoce el nivel real 
alcanzado por el agua al inicio de la falla, se 
deben definir el volumen descargado (V, ) por 
la brecha de falla, así como el tirante de agua 
sobre el fondo de la brecha H,. En la figura 1, H , 
es la altura del fondo de la brecha y H, la carga 
hidráulica sobre la corona. 

Otra información que es necesaria para 
ratificar las estimaciones realizadas es: (1) la 
longitud de cortina (L.), la distancia horizontal 
total en metros medida a lo largo del eje de la 
cortina en su cresta, es decir, entre los confines 
o bordes de su corona, que permite verificar el 
ancho máximo (b,,) de la brecha, la cual debe 
ser menor que L; y (2) el volumen de la cortina 
(V.), el cual debe ser mayor que la estimación 
del volumen erosionado (V,) de material en la 
brecha de falla. También es necesario conocer 
el tipo de cortina, para la aplicación de las 
ecuaciones empíricas más recientes. 


Recopilación exhaustiva y selección de 
ecuaciones empíricas 


Una búsqueda exhaustiva en todas las refe- 
rencias citadas condujo a la integración de una 
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Falla por desbordamiento 


Falla por tubificación 


Fondo del cauce 


PP ba —] 


Fondo del cauce 


Figura 1. Elementos geométricos generales de una brecha de falla (Froelich, 2008; Xu y Zhang, 2009). 


lista de 52 ecuaciones empíricas. Una primera 
selección consistió en eliminar las aplicables a 
presas naturales originadas por deslizamientos 
u obstrucción de ríos (Walder y O'Connor, 
1997); otra eliminación se basó en la similitud de 
fórmulas, en relación con sus coeficientes y/o 
exponentes, dejando sólo las representativas 
de un tipo de ecuación. A continuación, con 
base en una muestra de presas o embalses 
analizados en México (Campos-Aranda, 2013), 
se definieron magnitudes representativas de 
alturas de cortina y capacidades para realizar 
contrastes numéricos y eliminar las fórmulas 
con resultados dispersos o claramente erró- 
neos. Lo anterior condujo a seleccionar sólo 
28 ecuaciones empíricas. Tales ecuaciones se 
exponen en seguida en orden cronológico de 
aparición, agrupadas en las tres estimaciones 
siguientes: b, T, y Q,. 


Lista de variables utilizadas 


Para facilitar la comprensión rápida de cada 
una de las 28 fórmulas empíricas que serán 
expuestas, a continuación se presenta la lista de 
sus variables involucradas. 


m 


HM = 


variable de control (i varía de 1 a 5). 


= ancho máximo de la brecha de falla, en 


metros. 

ancho medio de la brecha de falla, en 
metros. 

ancho mínimo o del fondo de la brecha 
de falla, en metros. 

constante empírica, función del V, 
adimensional. 

tipo de cortina de corazónimpermeable. 


= tipo de falla por desbordamiento. 


aceleración de la gravedad igual a 9.81 
m/s?. 

altura de la brecha o profundidad desde 
la corona de la cortina, en metros. 
altura de la cortina, distancia desde el 
cauce hasta su corona, en metros. 

carga hidráulica sobre la corona de la 
cortina, en metros. 

altura del fondo de la brecha sobre el 
cauce, en metros. 

altura hidráulica, distancia del cauce 
a la elevación máxima del agua 
almacenada, en metros. 

tipo de cortina homogénea o de 
materiales graduados (mixta). 
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H,  = altura de cortina de referencia igual a 
15 metros. 

H,, = tirante de agua sobre el fondo de la 
brecha, en metros. 

k,, k, = factor por tipo de falla, adimensional. 

L.  = longitud de la cortina, medida a lo 
largo de su eje en su corona, en metros. 

PC = tipo de cortina de pantalla de concreto. 

Q, = gasto máximo de descarga, en m?/s. 


Op = gasto máximo de descarga, según curva 
envolvente, en m*/s. 

V  =mvolumen de almacenamiento O 
capacidad máxima de agua, con el 
sobre almacenamiento, en m?. 

V. = volumen de la cortina, en m?. 

V  = volumen de material erosionado de la 
cortina, en m. 

V,, = volumen descargado por la brecha de 
falla, en m?. 

v,  = variable de control (i varía de 1 a 5). 

Tb  = tipo de falla por tubificación. 

T,  = tiempo que dura la falla, en horas. 

T.  = tiempo de falla de referencia igual a 

una hora. 

taludes de los lados de la brecha de 

falla, adimensional. 


N 
Il 


Estimación del ancho de la brecha de falla 


Zagonjolli (2007), para el caso de 23 presas 
con falla por desbordamiento y cortinas de 
materiales térreos (homogéneas, mixtas y de 
enrocamiento) con o sin revestimiento en sus 
taludes y valores observados de b que variaron 
de 5 a 170 metros, encontró que la ecuación 
(2) siguiente tiene un 30% menos de error 
medio cuadrático (EMC) que la fórmula (1) 
originalmente propuesta por el US Bureau of 
Reclamation en 1988: 


b=3.0-H, (1) 

b=4.6:H, -5.7 (2) 

Por otra parte, Von Thun y Gillette (1990), 
citado en NWS (2002) y por Wahl (2004), 


analizaron 57 fallas de presas y propusieron la 
fórmula: 


b=2.5-H,+C, (3) 


en la cual C, es una constante función del 
volumen almacenado (V,) en millones de m* 
(Mm*), con los valores siguientes: 6.10 con V, 
< 1.233 Mm; 18.3 cuando V, fluctúa entre 1.233 
y 6.165 Mm; 42.7 cuando V, varía de 6.165 a 
12.330 Mn, y 54.9 con V, > 12.330 Mx*. 

Reed y Halgren (2011), basándose en el 
análisis de incertidumbre realizado por Wahl 
(2004) a las ecuaciones empíricas, recomiendan 
las encontradas por Froehlich (1995b), para b 
es: 


b=0.1803-k, -(H,) (4) 


e 
en donde k, es el llamado factor por tipo de 
falla, con un valor de 1.40 para desbordamiento 
y 1.00 para tubificación. H, es la altura o 
profundidad de la brecha en metros (ver figura 
1D). 

Froehlich (2008) estableció y presenta 
un registro de 74 fallas de presas, los cuales 
incluyen 69 datos de anchos medios de la 
brecha. Para tales valores encontró la expresión: 


b=0.270-k, + H00% .170320 6) 


en la cual b está en metros; k, es adimensional y 
vale 1.30 para fallas debidas a desbordamiento 
y 1.00 para los otros tipos, tubificación por 
ejemplo; Ves el volumen de agua almacenado 
arriba del fondo de la brecha en m? y que 
por lo tanto define el hidrograma de egresos. 
Froehlich (2008) también especifica que el talud 
z (figura 1) es aproximadamente 1.00 en fallas 
por desbordamiento y de 0.70 en las de otros 
tipos. 

Xu y Zhang (2009) presentan dos ecuaciones 
empíricas, una completa y otra simplificada, 
para las cinco características siguientes: H,, b, y b, 
Ty O, Como utilizan variables adimensionales, 
primero se debe estimar H,. En sus ecuaciones, 
H, = 15 metros es la altura de referencia de 
las cortinas, cuyos tipos pueden ser: (1) con 
corazón impermeable, (2) con pantalla de 
concreto, y (3) homogéneas o de materiales 
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graduados (mixtas); además establecen como 
variable predictiva principal la erosionabilidad 
de la cortina, clasificándola en alta, media y 
baja. En algunas ecuaciones se toma en cuenca 
el tipo de falla: desbordamiento o erosión por 
filtración (tubificación). Para la obtención de sus 
ecuaciones empíricas siguieron un proceso de 
selección de variables, basado en el coeficiente 
de determinación ajustado. Sus primeras tres 
ecuaciones completas y simplificadas permiten 
estimar H,, b,, y b y corresponden a las números 
(6), (7) y (8) para las primeras, y a las (9), (10) 
y (11) para las segundas, todas mostradas en el 
cuadro 1. 


Estimación del tiempo de formación de la 
brecha de falla 


MacDonald y Langridge-Monopolis (1984) son 
los primeros en proponer una relación empírica 
gráfica en papel logarítmico para estimar T, 
en horas, en función del volumen de material 
erosionado (V,) de la cortina en m'. Además, 
definen el llamado factor de formación de la brecha 
como el producto de la altura hidráulica por el 
volumen de agua descargado (H,'V, ), el cual 
permite estimar V. Wahl (2004) ha expuesto 
estas relaciones como ecuaciones potenciales, 
éstas son: 


T; =0.0179 ad (12) 


en la cual V, se estima con las expresiones 
siguientes para cortinas de tierra homogéneas 
y de materiales graduados (mixtas) y de otros 
tipos (enrocamiento y terraplenes con pantalla 
de concreto o con corazón impermeable): 


(V.),, =0.0261-(H, A il 


(13) 


(14) 


(Wales = 0.00348-(H, V, ld 


D 


Nuevamente, Reed y Halgren (2011) 
adoptan la fórmula propuesta por Froehlich 
(1995b), según resultados de Wahl (2004): 
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T, =0.00254-(H,)%-( 


V, 


0.53 
w ) 


(15) 


Froehlich (2008) encuentra con análisis 
dimensional y procesamiento de 23 datos la 
fórmula siguiente: 


1/2 
A 


2 
b 


T, =0.0176- (16) 


Por último, las ecuaciones completa y 
simplificada que encontraron Xu y Zhang 
(2009) para el tiempo de falla son las números 
(17) y (18) del cuadro 1, cuyo T, es de una hora. 


Estimación del gasto máximo de descarga 


Hagen (1982), citado por Wahl (2004) y 
Zagonjolli (2007), analizó 18 fallas de presas por 
desbordamiento y encontró que el producto 
de H, por V, o V, , llamado factor de presa, se 
relaciona con el gasto máximo de descarga Q,. 
Los gastos máximos observados variaron de 
730 a 83 000 m*/s, mientras que el factor de 
presa fluctuó de 2.90 a 48 000 Mm*. Su ecuación 
envolvente propuesta es: 


Qp=0.5404-(H,-V,,) 


w (19) 

Walder y O'Connor (1997) señalan que las 
llamadas curvas envolventes definen el límite 
superior del Q, en las gráficas, que muestran 
su valor en función de H,, V,, o de (H,'V,,) y 
tienen como propósito único permitir hacer 
una estimación rápida y conservadora del Qp. 
Las curvas envolventes no tienen significado 
físico y sólo han sido establecidas para mostrar 
los límites superiores de las observaciones 
históricas. 

MacDonald y Langridge-Monopolis (1984) 
analizaron la falla de 42 presas, 30 con cortinas 
homogéneas, y 12 con cortinas mixtas y de 
enrocamiento con pantallas de concreto. La 
altura de ellas varió de 6 a 93 metros. Definieron 
dos ecuaciones de regresión entre Q, y el factor 
de formación de la brecha; la segunda es la 
curva envolvente: 
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Q, =1.154-(H, *V, 


0.412 
A 


(20) 


Qp=3.850-(H, -V,,) 


, (21) 

Costa (1985), también citado por Wahl 
(2004) y Zagonjolli (2007), analizó 31 fallas 
históricas de presas, cuya altura varió de 18 a 
83.8 metros y volumen de almacenamiento al 
tiempo de la falla, fluctuando de 3 800 m* a 700 
Mrm?. No se hizo distinción en el modo de falla 
ni en el tipo de cortina. Sus dos ecuaciones de 
menor incertidumbre según el estudio de Wahl 
(2004) son: 


(22) 


w 


Qp=2.634-(H,V,) 


pp (23) 


Q, =0.981-(H,, *V,, 

Zagonjolli (2007) encontró gráficamente 
que no existen diferencias significativas entre 
las ecuaciones (20) y (23). Froehlich (1995a,b) 
estudió 22 fallas históricas de presas, cuyas 
alturas fluctuaron de 3.7 a 86.9 metros; sus 
volúmenes V,, variaron de 92 500 m*a 660 Mim*. 
La fórmula encontrada es: 


Q, =0.607-(H,,) "-(V, 


0.295 
w) 


(24) 


Froehlich (1995a) cita como referencia 
histórica la fórmula propuesta por el US Soil 
Conservation Service en 1985 para cortinas de 
materiales térreos y valores de H, mayores de 
31.4 metros: 

Q, =16.60-(H,)* (25) 

Walder y O'Connor (1997) tratan el tema de 
las presas naturales, es decir, las formadas por 
los procesos que originan un almacenamiento 
de agua al obstruir un cauce; los más comunes 
son los deslizamientos de ladera y los flujos 
con escombros (tallos y ramas). Estos autores 
presentan tres ecuaciones de regresión y 


sus respectivas curvas envolventes para 
presas originadas por obstrucción y para 
las construidas ex professo. Para las presas 


construidas sus fórmulas son: 


Q, =2.50-(H, (26) 
Q, =1.16-(V,,) (27) 
Q, =0.61-(H, AN (28) 
Qp=2.90-(H, VW, (29) 


Walder y O'Connor (1997) indican que 
sus fórmulas están basadas en los datos 
recolectados por Costa (1988). Estos autores 
señalan que probablemente las estimaciones 
del Q, obtenidas con las ecuaciones de regre- 
sión sean menores que el gasto real descarga- 
do por la brecha de falla debido a que muchas 
mediciones se realizaron aguas abajo de la 
presa y tales valores no fueron corregidos 
por traslación y atenuación de la onda de la 
creciente. Esto último también justifica el uso 
de las curvas envolventes para acotar el valor 
del Q,. Las últimas dos ecuaciones de Xu y 
Zhang (2009) completa y simplificada son las 
números (30) y (31) del cuadro 1 y permiten la 
estimación de O, 


Factores de mayor influencia en las 
estimaciones 


En la mayoría de las ecuaciones simplificadas 
encontradas por Xu y Zhang (2009), el factor 
más importante que influye en la estimación 
es la erosionabilidad de la cortina. Una 
erosionabilidad alta permite un desarrollo 
rápido de la brecha de falla, lo cual conduce 
a dimensiones mayores de ésta, a un gasto 
máximo de descarga más elevado y un 
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tiempo de formación más breve. También 
se Observa que el llamado cociente de forma 
del almacenamiento me H,,) tiene un papel 
preponderante en todas sus ecuaciones, 
excepto en las dos primeras, que estiman la 
profundidad de la brecha (H,). 

En relación con el cociente de forma del 
almacenamiento existen dos tipos de presas o 
embalses: la llamadas bajas con altura reducida 
de cortina y gran almacenamiento, es decir, 
embalses de planicie, que conduce a valores 
de vi A H,, grandes; y las presas altas con una 
altura considerable de cortina y volumen de 
almacenamiento reducido, esto es, embalses 
encañonados, que originan magnitudes de 
Jal el H., pequeñas. Entonces, para dos presas 
con el mismo volumen de almacenamiento, 
pero una alta y otra baja, la primera generará 
un gasto de descarga mayor debido a la mayor 
energía potencial, pues como se observa en 
las ecuaciones (30) y (31), el gasto Q, decrece 
conforme Vd / H,, se incrementa. T, decrece 
con el incremento de Q, y aumenta conforme 
vi 3 /H,, se acrecienta (ver ecuaciones (17) y 
(18). 

Varios investigadores han encontrado que 
el gasto máximo puede ocurrir antes de que 
el ancho de la brecha de falla sea máximo. Xu 
y Zhang (2009) indican que en presas altas es 
muy probable que el gasto máximo de descarga 
ocurra antes que en las presas bajas y que en 
estas últimas el ancho de la brecha de falla será 
mayor debido a que la erosión lateral continúa 
con un nivel en el embalse relativamente alto. 

También se ha encontrado que las fallas por 
desbordamiento originan anchos de brecha 
mayores que los rompimientos por erosión 
debida a filtraciones (tubificación); es probable 
que lo anterior esté asociado principalmente a 
que las fallas por desbordamiento se originan 
por lo común por crecientes que aportan 
grandes volúmenes de agua al embalse, los 
cuales elevan su nivel más allá de la corona de 
la cortina, aumentando su cociente Y > Hu 
por el contrario, en las fallas por tubificación 
se libera agua del embalse y los niveles en el 


almacenamiento son menores tanto en el inicio 
del proceso como después del colapso del túnel 
de erosión para formar la brecha de falla. 


Intervalos de incertidumbre asociados con 
varias ecuaciones empíricas 


En el cuadro 2 se presentan las estimaciones de 
Wahl (2004) relativas a la incertidumbre de 11 
de las 28 ecuaciones empíricas seleccionadas. 
Se indica el error medio de la estimación en 
ciclos logarítmicos (2), el ancho de la banda 
de incertidumbre (+25), también en ciclos 
logarítmicos, siendo S, la desviación estándar 
de los errores y el intervalo de variación de 
la estimación o intervalo de confianza (1 
de C.), los cuales corresponden a 10% y a 
107*2., mismos que al ser multiplicados por 
el valor obtenido con la ecuación empírica 
respectiva definen aproximadamente la banda 
de confianza al 95%. Para las 10 ecuaciones 
empíricas propuestas por Xu y Zhang (2009) 
se exponen en el cuadro 2, en sus dos últimas 
columnas, sólo los intervalos de variación 
establecidos por tales autores. 


Agrupamiento de las ecuaciones empíricas 


Se integraron tres grupos con las 28 ecuaciones 
empíricas que son aplicables a presas cons- 
truidas, por requerir sólo información de altura 
del nivel del agua en el almacenamiento y 
volumen del mismo antes de la falla. El primer 
grupo lo forman siete ecuaciones seleccionadas 
que no fueron contrastadas por Wahl (2004) 
y que por lo tanto no tienen intervalos de 
confianza para sus estimaciones; el segundo 
grupo de 11 ecuaciones son las reportadas 
en el cuadro 2 y que fueron estudiadas por 
Wahl (2004). El último grupo lo constituyen 
las 10 ecuaciones empíricas encontradas por 
Xu y Zhang (2000), consideradas como las 
más confiables debido a la información que 
incorporan en las variables de control (B, y C, 
del cuadro 1) sobre la erosionabilidad de la 
cortina y su tipo, así como por su forma de 
romperse o fallar. 
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Cuadro 2. Estimaciones de la incertidumbre de las ecuaciones empíricas indicadas, 
según Wahl (2004), y Xu y Zhang (2009). 


Número de ecuación Intervalo de Número de rcuación Intervalo de 


(variable por estimar) á sea confianza al 95% (variable por estimar) confianza al 95% 
(4) (b) +0.01 + 0.39 0.40-2.40 (6) y (9) (H, < 50 m) 0.79-1.20 
(3) (b) +0.09 = 0.35 0.37-1.80 (6) y (9) (H, > 50 m) 0.74-1.25 
(12) (T) -0.21 + 0.83 0.24-11.0 (7) (0,1) 0.44-2.27 


=0:22 + 0.64 0.38-7.30 (10) (b,,) 0.45-2.22 


(Q4) (Q,) 0.04 + 0.32 0.53-2.30 (8) (b) 0.43-2.36 
(23) (Q,) 0.05 + 0.72 0.17-4.70 (11) (b) 0.44-2.30 
(25) (Q,) +0.,13 + 0.50 0.23-2.40 (17) (Tr) 0.30-3.35 
(20) (Q,) +0.,13 + 0.70 0.15-3.70 (18) (T) 0.34-2.97 
(19) (Qp) +0,43 + 0.75 0.07-2.10 (32) (Q,) 0.30-3.35 
(21) (Qp) +0.64 + 0.70 0.05-1.10 (33) (Q,) 0.33-3.02 
(22) (Qp) +0.64 + 0.72 0.04-1.22 = = 
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Aplicaciones numéricas 


Casos hipotéticos establecidos y sus 
resultados 


Para ilustrar los resultados de las 18 ecuaciones 
empíricas que son aplicables cuando sólo existe 
la información básica que se ha citado (H, y V.), 
se han formulado dos aplicaciones numéricas 
hipotéticas. La primera corresponde a una 
presa relativamente pequeña, con altura de 
cortina (H,) de 13 metros, altura del agua a nivel 
de conservación (NAN) de 10 metros y cuya 
falla por desbordamiento originó un tirante 
de agua sobre la corona de 1 metro, entonces 
H,, es de 14 metros, ya que se supone que la 
brecha de falla llegó al desplante de la cortina, 
en el nivel del cauce. El volumen almacenado 
a nivel de conservación (V,) es de 15 Mn y al 
inicio de la falla (V,,) de 30 Mm”. La cortina es 
homogénea de material arcilloso, que puede 
ser considerado fácilmente erosionable. 

La segunda aplicación es la relativa a una 
falla por tubificación en una presa de mediana 
altura, cuyos datos son: H, = 45 m, H,, = 35 
m, altura del NAN = 40 m, V, = 140 Mm* y 
V,, = 80 Mrw*. Su cortina tiene corazón imper- 
meable y puede ser considerada resistente a la 
erosión (erosionabilidad media). 


En la primera aplicación se emplean k, =1.40 
y k, = 1.30, en las ecuaciones (4) y (5); el tipo de 
cortina es H para estimar el volumen erosio- 
nado con la expresión (13) y se usan las ex- 
presiones (12), (15) y (16) para estimar el T, Se 
proponen dos análisis: el primero utilizando 
los datos fácilmente disponibles de la presa, es 
decir, altura del NAN y volumen de conserva- 
ción, por lo cual H, = 10 metros y V,, =15 Mir?. 
En el segundo análisis se emplea el nivel antes 
de la falla y el volumen almacenado respectivo, 
entonces H, = H,,=14 metros y V,,=30 Mim*. Los 
resultados de la aplicación de las 18 ecuaciones 
empíricas en cada análisis se muestran en el 
cuadro 3, así como los valores adoptados en 
cada uno. 

Para la segunda aplicación se tiene que 
k, = k, = 1.00; en las ecuaciones (4) y (5), el 
tipo de cortina es OT para estimar el volu- 
men erosionado con la expresión (14) y se 
usan las expresiones (12), (15) y (16) para 
estimar el T, De nuevo se proponen dos 
análisis: el primero utilizando los datos 
fácilmente disponibles de la presa, es decir, 
altura del NAN y volumen de conservación, 
por lo cual H, = 40 metros y V,, = 140 Mm; 
en el segundo análisis se emplea el nivel antes 
de la falla y el volumen almacenado respectivo, 
entonces H, = H,,=35 metros y V,, = 80 Mm*. Los 
resultados de la aplicación de las 18 ecuaciones 
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Cuadro 3. Resultados de las 18 ecuaciones empíricas indicadas, en la primera aplicación numérica: 


falla por desbordamiento y cortina homogénea erosionable. 


Primer análisis: 
H, = 10.0 m y V,, =15.0 Mm? 


Segundo análisis: 
H, = 14.0 m y V, = 30.0 Mm? 


Número de Magnitud Intervalo de Número de Magnitud Intervalo de 
ecuación estimada confianza ecuación: estimada confianza 
(4) b=77.4 30.9-185.6 (4) b= 102.9 41.2-247.0 
(3) b=79.9 29.6-143.8 (3) b=89.9 33.3-161.8 
(2) b=40.3 (2) b =58.7 
(5) b=76.1 (5) b=96.3 
Bao 750 m Dadop = 95.0 m 
(13) V, = 50 589 (13) V,= 111 665 
(12) US 0.22-10.15 (12) T,=1.23 0.30-13.54 
(15) T,= 2.03 0.77-14.84 (15) T,=2.17 0.82-15.83 
(16) Ia, (16) T,=2.20 
Mondo = 1-75 h (Thadop = 2-00 h 
(4) Q, =1381 732-3 176 (24) Q,=2571 1 363-5 914 
(23) Q, =2 665 453-12 524 (23) Q, = 4 106 698-19 299 
(25) O), = 1105 270-2 821 (25) Q,=2190 504-5 256 
(20) Q, =2 696 404-9 976 (20) Q, = 4121 618-15 247 
(19) Qp =6 619 463-13 899 (19) Qp = 11 075 775-23 257 
Q1) Qp = 8827 441-9710 (Q1) Qp = 13 478 674-14 826 
(Q2) Qp =10 425 417-12 719 (22) Qp =16 400 656-20 008 
(26) Q, = 547 (26) Q,=1 202 
(27) Q, =2320 (27) Q,=3191 
(28) Q, =2 000 (28) Q,=3114 
(29) Qp =9 509 (29) Qp = 14 805 
(Qaciop =2 000 m*/s 425-6 500 (QOaio, =3 500 m/s 650-10 000 
(Qp), cop = 10 000 m?/s 450-12 000 (Qp), dop = 15 000 n?/s 700-20 000 


empíricas en cada análisis se muestran en el 
cuadro 5, así como los valores adoptados en 
cada uno. 

En los cuadros 4 y 6 se exponen los resulta- 
dos de la aplicación de las ecuaciones comple- 
tas y simplificadas de Xu y Zhang (2009) para 
la primera y segunda aplicaciones numéricas. 


Análisis de los resultados 


En relación con los resultados concentrados 
en los cuadros 3 a 6, correspondientes a las 
dos aplicaciones numéricas planteadas y sus 
respectivos enfoques de uso de la información 
disponible, es destacable que todos ellos 


muestran una similitud en órdenes de mag- 
nitud, en particular las estimaciones del ancho 
de la brecha y del gasto máximo de descarga. 
Considerando que las estimaciones mostradas 
en los cuadros 4 y 6 son las más confiables, es 
importante citar que el uso de las 18 ecuacio- 
nes empíricas seleccionadas conduce a valores 
de b y T,aproximados en la primera aplicación, 
no así con su O; lo contrario ocurre en la 
segunda aplicación numérica. 

Respecto a los resultados de las ecuaciones 
completa y simplificada de Xu y Zhang 
(2009), mostrados en los cuadros 4 y 6 para 
cada aplicación numérica, se observa para la 
primera el efecto de la erosionabilidad alta, la 
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Cuadro 4. Estimación de las magnitudes asociadas con el rompimiento de una presa con cortina homogénea de tierra, de alta 
erosionabilidad y falla por desbordamiento, con las ecuaciones de Xu y Zhang (2009). 


Datos generales: H, =13.0 m; H,,= 14 m; V, = 30.0 Mm* 


Ecuación completa Ecuación simplificada 
Número de Magnitud Intervalo de Número de Magnitud Intervalo de 
ecuación estimada confianza ecuación estimada confianza 
(6) H,=135m 10.6-16.2 (9) H,= 187 m 10.8-16.4 
(7) b,,=126.7 m 55.8-287.7 (10) b,, 144.7 m 65.1-321.2 
(8) b=97.1m 41.8-229.2 (11) b=121.1m 53.3-278.6 
(17) 1 M2 h 0.52-5.76 (18) T,=1.65 h 0.56-4.90 
(30) Q, = 4 483 n?/s 1345-15 018 (31) Q, = 5653 m/s 1865-17 071 


Cuadro 5. Resultados de las 18 ecuaciones empíricas indicadas en la segunda aplicación numérica: falla por tubificación 


y cortina resistente a la erosión. 


Primer análisis: Segundo análisis: 
H, = 40.0 m y V,, =140.0 Mm? H, = 35.0 m y V,, = 80.0 Mm? 
Número de Magnitud Intervalo de Número de Magnitud Intervalo de 
ecuación estimada confianza ecuación estimada confianza 
(4) b=146.9 58.8-352.7 (4) b=119.8 47.9-287.5 
(3) b=154.9 57.3-278.8 (3) b=142.4 52.7-256.3 
(2) o = 1743 Q) (0 = 1535 
(5) b=126.5 (5) b=105.2 
BA =150.0 m Da = 130.0 m 
(14) V, =703 130 (14) V, = 389 545 
(12) T,=2.41 0.58-26.46 (12) SA 0.47-21.34 
(15) T,=191 0.72-13.92 (15) T,= 1.60 0.61-11.67 
(16) T,= 1.66 (16) TILA 
To, dop =2-00 h To, dop = 1.65 h 
(Q4) Q, = 14 889 7891-34 244 (Q4) Q,=10 697 5669-24 602 
(Q3) Q, = 12 187 2072-57 280 (23) Q, =9 109 1 549-42 813 
(25) Q, =15 273 3513-36 655 (25) Q,= 11 90 2 744-28 632 
(0) Q,= 11 980 1797-44 326 (0) Q, = 9004 1351-33 314 
E (19) Qp =40 440 2 831-84 924 (19) Qp =28 595 2 002-60 050 
a Q1) Qp =39 081 1 954-42 989 (1) Qp =29 393 1 470-32 332 
El (Q2) Qp =51 265 2051-62 543 (Q2) Qp =37 789 1512-46 103 
a (26) Q, = 14 020 (26) Q, = 10 258 
5 (27) Q, = 6 482 (Q7) Q, =5011 
E (Q8) Q, = 9 485 (28) Q, =7041 
: (29) Qp = 45 094 (29) Qp =33 472 
= (0), dop 13 500 m/s 2 500-40 000 (Q), dop — 10 000 m/s 2 000-30 000 
(Qp), dop — 45 000 m?/s 2 250-62 000 (Qp), dop — 35 000 m/s 1500-45 000 
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Cuadro 6. Estimación de las magnitudes asociadas con el rompimiento de una presa con cortina de enrocamiento con corazón 


impermeable, de erosionabilidad media y falla por tubificación, con las ecuaciones de Xu y Zhang (2009). 


Datos generales: H, = 45.0 m; H,,= 35 m; V,, = 80.0 Mm? 


Ecuación completa 


Ecuación simplificada 


Número de Magnitud Intervalo de Número de Magnitud Intervalo de 
ecuación estimada confianza ecuación estimada confianza 
(6) HA,=41.7m 33.0-50.1 (9) H,=41.0m 32.4-49.2 
(7) b,,=138.0 m 60.7-313.2 (10) b,,=116.3 m 52.3-258.1 
(8) b=110.4 m 47.5-260.5 (11) b=86.6 m 38.1-199.1 
(17) T,=3.21 h 0.96-10.75 (18) T,= 3.09 h 1.05-9.18 

(30) O, =1579 m?/s 4 738-52 909 (31) O, =11571 m*/s 3 819-34 946 

cual conduce a estimaciones mayores en b y Tercera: en el proceso de selección de 


Q,, Lo inverso ocurre en la segunda aplicación 
numérica, quizás debido a la erosionabilidad 
media. 


Conclusiones 


Primera: la selección expuesta de ecuaciones 
empíricas que permiten estimar las magni- 
tudes asociadas con la falla de una presa por 
rompimiento de su cortina de materiales 
térreos, es decir, el ancho medio de brecha (b) 
en metros, el tiempo de formación de ésta (1) 
en horas y el gasto máximo descargado (Q,) en 
m?/s, constituye una primera aproximación en 
los estudios de seguridad de presas o embalses, 
encaminados a definir los planes de emergencia 
y las estrategias de mitigación de riesgos en las 
áreas localizadas aguas abajo de éstas. 

Segunda: a juzgar por los resultados con- 
centrados en los cuadros 3 a 6, para las dos 
aplicaciones numéricas formuladas, éstos per- 
miten establecer los valores probables de las 
magnitudes b, T, y Q,, así como sus intervalos 
de variación; con base en tales estimaciones, 
se pueden plantear los análisis de sensibilidad 
que deberán ser abordados con los modelos 
hidráulicos de estimación del hidrograma de 
egresos y su tránsito respectivo hacia aguas 
abajo, para definir gastos, niveles y tiempo 
de llegada de la inundación que se estima, 
originada por rompimiento o falla de una presa 
con cortina de tierra o de enrocamiento. 


resultados se deben tener en cuenta las dos 
siguientes recomendaciones: (1) se considera 
que las 10 ecuaciones empíricas desarrolladas 
por Xu y Zhang (2009), mostradas en el cuadro 
1, son las más confiables, debido a que fueron 
generadas con un proceso de selección de 
regresores e incorporan en sus variables de 
control (vw) información relevante sobre la 
cortina, su tipo de falla y la erosionabilidad 
estimada; (2) del resto de las 18 ecuaciones 
empíricas, tienen mayor confiablidad las 11 
que presentan intervalos de variación, pues 
fueron las que mostraron menor incertidumbre 
de acuerdo con el trabajo de Wahl (2004). 
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